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Abstract 

Within the MSSM we calculate the electroweak 1-loop box contribu- 
tions to the processes e + e~ — > Z°h° and e + e~ — > A°h°. We present 
detailed results for c. m. energies y/s = 200 GeV and y/s = 500 
GeV as well as for tan (3 = 2 and tan /3 = 50. The box contributions 
to the process e + e~ — > are, depending on y/s and tan/3, of 

the order -2 to -20% and to e+e" -> A°h° of 2 to 10%. 



1 Introduction 



Experiments at the forthcoming upgrade of the LEP collider to LEP2 will continue the search for 
Higgs bosons in the mass range up to 100-110 GeV depending on the energy and luminosity 
at this machine. Besides the search for the Standard Model Higgs boson, numerous studies are 
devoted to the Higgs bosons of the minimal supersymmetric standard model (MSSM) g, the 
most predictive framework beyond the Standard Model, where at least one light neutral scalar 
with mass ~ 130 GeV is expected. This mass range can be fully explored at new e + e~ colliders 
with energies up to 500 GeV 0, [/]. Detailed studies of Higgs production and decay processes 
are required to detect the possible signals of Higgs bosons and to distinguish as far as possible 
the origin of a produced scalar particle. The most promising process to search for the standard 
Higgs boson is the bremsstrahlung type reaction e + e~ — > ZH at energies > Mh + Mz- For 
Higgs bosons of the MSSM the basic production channels are 

e+e" -> Zh° (Z) 

and 

e + e~ -> A°h° (A) 

where h° is the lightest CP even and A the CP odd neutral particle of the MSSM Higgs 
spectrum. 

For an accurate discussion of the production cross sections one has to include radiative 
corrections. For the standard Higgs complete 1-loop calculations for e + e~ — ► ZH are available 
H. For the corresponding processes in the MSSM no complete 1-loop calculation has been 
performed so far. Besides the approximate methods of the effective potential approach M and 
the renormalization group treatment [H it is of particular interest to also have a complete 
diagrammatic calculation: 

- it allows for all virtual particle effects within the MSSM without restrictions to masses 
and mixings, 

- it takes into account all momentum dependent effects in 2- and 3- point functions, 

- it provides a reference frame for checking the quality of simpler compact approximations 
which might be useful for practical applications. 

The complete set of 1-loop diagrams contributing to the MSSM neutral Higgs production pro- 
cesses is very extensive, and until now only one calculation has been carried out || based on the 
following building blocks: the set of self energy corrections to the vector boson and Higgs prop- 
agators and wave functions, the vertex corrections to the e + e~~Z, ZZh, ZAh vertices, and the 
corresponding counter terms in the on-shell renormalization scheme. The box diagram contri- 
butions have not yet been considered so far, and they are also not part of the phenomenological 
studies for LEP2 and higher energy colliders [0. |J. 
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In this paper we provide the box contributions to both e + e~ — > Zh and e + e~ — > A/i. Box 
diagrams are required for theoretical reasons to render the amplitudes gauge independent, and 
for practical reasons since they are in general of non-negligible size. The computation in || on 
the basis of 2- and 3-point functions was performed in the 't Hooft-Feynman gauge with £ = 1. 
In this specific gauge, the box contributions are finite and independent of the details of the 
renormalization. Hence, we proceed to derive the box contribution in the R^=i gauge as well, 
thus providing the missing ingredients as an independent, finite and scheme independent block 
for completing the amplitudes and cross sections. 



2 Cross Sections 

The Born diagrams contributing to process (Z) and (A) are displayed in Fig. 1 together with 
the momentum assignment. p 3 always denotes the outgoing h° momentum, whereas p^ is the 
momentum of the outgoing Z or A , respectively. e z denotes the polarization vector of the Z . 
Furthermore, we use the invariant variables 

s= (pi+ P2) 2 , t = (pi - p 3 ) 2 , u={px- p 4 ) 2 

and the scattering angle 9 = ^(e~, h°) in the c. m. system. The masses of e + and e~ have been 
neglected for the entire calculation. 

Including the box diagrams (Fig. 2 and 3) one has to take into account their interference 
with the corresponding Born amplitudes. The differential cross section in the c. m. system can 
be written in the following way: 



da z J\(s,M 2 ,M 2 ) 



dcos# 327rs 5 



da A J\(s,M 2 A ,M 2 ) 



d cos 9 327rs 5 



n+j:T i A ) (2) 



with A(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz). The index n counts the diagrams of Fig. 2, 
and / those of Fig. 3. Thereby it is understood that the summation is performed also over the 
various x° an d X + states as well as over el and e#. 

T z A denote the lowest order contributions following from the Born amplitudes 

M% A = v( P2 ) 0% ufa) (3) 

with 



q 2 My sin(/? — a) r. e . „ 

C W {S +vl z ) 



W (4) 



for process (Z), and 

\(9v + 9a)Vl + (g e v - gl)V R ] (j9 3 - p,T (5) 



o g 2 cos((3 - a) 

2 4( S -M|) 7 ^ 
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for process (A), where 

9v = -\ + *l, 9a = -\ (6) 
are the standard NC couplings, g = e/s w , and a, [3 the mixing angles of the Higgs sector. 

Averaging over the initial helicities for unpolarized e ± and summing over the Z polarizations 
a in case of (Z) yield the following expressions 

T z = \T,^(^OUiK) (7) 

n = ^(j/ 2 o^i0°j (8) 

which result in the explicit forms of the Born cross sections: 



AO% - 1 ^ ^ - a) ^ ^ M t M } ' fa 2 + <&') 2 (« » + 2 S Mf - MfM*) (9) 



dcosfl 128tts 2 C 4 ^ 7 (s-Mf) 2 



The interference terms in ([[]), (0) can be expressed in a way analogous to (|7|.|8|): 

Tz = \ ^Re\Tr(^O z ^ z )] (11) 



T\ = \ Re [Tr (j/ 2 O ^ 



:i2i 



with the operators corresponding to the diagrams of Fig. 2 and 0^ to those of Fig. 3. 

For the discussion of the various contributions it is useful to distinguish between the genuine 
SUSY diagrams with virtual selectrons (se) and sneutrinos (snu) and the residual diagrams of 
the standard model with 2 Higgs doublets (2hdm): 
Figure 2: 

Of = diagrams(l, 2, 3) + (h Z) 

C| nu = diagrams(4, 5, 6) + (h <-> Z) (13) 
e>| hdm = diagrams(7, 8) + (/i «-> Z) + diagrams(9, 10) 



Figure 3: 



= diagram 2 + (h ^ A) + diagram 4 
O s ™ = diagram 3 + (h «-> A) + diagram 5 (14) 
C» 2hdm = diagram 1 + (/i «-> A) 



including summation over the x°, X + ; an d e L « states. 
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Our conventions for the Feynman rules are adopted from ref. and for fermion number 
violating processes from ref. (KI . The calculations were done using the algebraic program FORM 



TTj] and the package FF JT/2J to evaluate the tensor integrals. 

According to the decomposition ( J13| , pH]) we can write for the cross sections: 

da? a da ? 



'Z,A 

dcos# 



' Z,A 



dcos# 



1 + 6 hox 



Z,A 



with the relative box contribution 



srbox r2hdm 

>Z,A — °Z.A "T "Z,A U Z,A- 



+ 07, A+0? 



In a similar way we express the integrated cross sections as 



with 



A 



OX4 



box 



a 



Z,A 



1 



A box 
^Z,A 



A2hdm 1 a se 



1 a snu 
+ ^Z,A 



(15) 



(16) 



(17) 



l Z,A — ^Z,A ~ ^Z,A 

where o\ A denote the Born terms following from (|9|) and ( |T0"D by integration over cos#. 



3 Results and discussion 

In this section we present and discuss the numerical influence of the box diagrams on the pro- 
duction cross sections. As concrete examples we have chosen two sets of Higgs masses where 
both Zh° and A°h° can be produced at LEP2 energies. Taking Ma = 90 GeV, one has = 40 
GeV for tan/? = 2 and M h = 89 GeV for tan/3 = 50. These mass values correspond to the 
lowest order MSSM Higgs spectrum. At this stage we have not yet included higher order effects 
since that would include a bunch of additional parameters like m t , scalar top masses and mixing 
parameters which are not directly related to the specific problems addressed in this note. The 
input required for calculating the box contributions is Mh and the mixing angle a, besides M^, 
tan/3, the slepton and chargino/neutralino parameters. In a full calculation M^ and a will be 
derived from the dressed (h°, H°) propagator system and thus are provided automatically in a 
comprehensive one-loop program. The simplified approach corresponds roughly to a situation 
where the scalar top system has masses around m t , i.e. -m^ = mf, when no large corrections 
to neutral Higgs masses and couplings occur. 

We discuss the two processes (Z) and (A) seperately: 
e+e- -»■ Zh°: 

In order to exhibit the size of the box contribution we show in Fig. 4 the relative correction 
<5| ox , eq. flTB]), in the differential cross section as function of the scattering angle 9 for two en- 
ergies, \/i = 200 GeV and 500 GeV. For the scalar leptons a universal mass mj = 100 GeV 
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and no mixing is assumed, and the x , X + system is fixed by choosing M = fi = 150 GeV. This 
corresponds to a mass of 76 GeV for the lightest x + and 44 GeV to the lightest x°- The figure 
shows that the contribution reaches already —5% at 200 GeV for large tan (3, and increases 
remarkably to —(20 4- 25)% at 500 GeV. It can be seen that the dominating piece is from the 
gauge and Higgs sector (2hdm), whereas the genuine SUSY boxes are smaller than 5%, also 
when we decrease the slepton masses to lower values close to the present experimental limit. 
This dominance of the gauge-Higgs part persists also at higher energies, as displayed in Figure 
5 showing the integrated cross section as function of the c. m. energy. The dependence on the 
SUSY breaking parameters M, n is visualized in Fig. 6, again for the integrated cross section: 
for increasing M and /x the genuine SUSY boxes become very small, as expected from decou- 
pling properties, and we are left with the gauge-Higgs subclass of diagrams. As a consequence, 
the box contribution to e + e~ — > Zh is not very sensitive to the details of the SUSY sector, and 
is qualitatively very similar to the Standard Model situation ||. 

e+e- -»■ A°h°: 

In Fig. 7 we show the relative correction 5^ ox (eq. [T^), with the same parameter choices as in 
Fig. 4. The contribution of the box diagrams to the process (A) are somewhat smaller than 
to the process (Z): at y/s = 200 GeV the relative correction is about 2% and reaches nearly 
10% at 500 GeV. Again the contributions of the gauge and Higgs sector dominate. The effect of 
diagrams containing genuine SUSY particles is less than 0.5-=- 1% for mr— 100 GeV and slightly 
larger for mj = 50 GeV. The dependence of the integrated cross section on the c. m. energy 
\/s is plotted in Fig. 8. As can be seen, the effects of genuine SUSY diagrams decrease with 
y/s. We also show the dependence on M and \i in Fig. 9: the behaviour is similar to the case 
of the Zh° production, i.e. the SUSY boxes become very small with increasing SUSY breaking 
parameters. 

In summary, the box diagrams yield a sizable contribution to the differential and integrated 
Higgs production cross sections, in particular in the energy domain of a next linear e + e~ collider. 
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Figures 

• Fig.l Born diagrams contributing to e + e~ — > Z°h° and e + e~ — > A°h°. 

• Fig. 2 Box diagrams contributing to e + e~ — > 

• Fig. 3 Box diagrams contributing to e + e~ — > A°h°. 

• Fig.4 Relative box contribution 5 b z ox for e+e" -> at >/s = 200 GeV, 500 GeV and 
tan/5 = 2, 50. 

• Fig. 5 Relative box contribution A b £ x for M = /x = 150 GeV. 
If curves are marked by '-se' etc. — Af A is plotted. 

• Fig.6 A b £ x &ty/s = 200 GeV, 500 GeV and tan/3 = 2, 50. 

• Fig. 7 Relative box contribution 6% x for e + e~ ^ A°h° at ^ = 200 GeV, 500 GeV and 
tan/3 = 2, 50. 

• Fig.8 Relative box contribution A b £ x for M = /i = 150 GeV. 

• Fig.9 A b X x at v/s = 200 GeV, 500 GeV and tan/3 = 2, 50. 
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Figure 2 
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Figure 5 
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Figure 7 
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eV -> A°h°, M susy =/x=1 50GeV, m A =90GeV, 1 OOGeV sf ermions 
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Figure 9 
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